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Abstract. We study the relationship between the H2 and CO abun- 
dances in simulated molecular clouds using a fully dynamical model 
of magnetized turbulence coupled to a detailed chemical network. We 

' find that the CO-to-Hb conversion factor for a given molecular cloud, 

the so-called X-factor, is determined primarily by the mean extinction 

■ of the cloud, rather than by its metallicity. Our results explain the 

discrepancy observed in low metallicity systems between cloud masses 
Q . derived from CO observations and other techniques such as infrared 

emission, and predict that CO-bright clouds in low metallicity systems 
should be systematically larger and/or denser than Milky Way clouds. 

7— ( 1 Introduction 

> , 

' Observed star formation takes place within giant molecular clouds (GMCs), so un- 
derstanding how these clouds form and evolve is a key step towards understanding 
star formation. The main chemical constituent of any GMC is molecular hydrogen 
(H2), but this is very difficult to observe in situ. For this reason it is common to 
use emission from carbon monoxide (CO) as a proxy for H2. 

In order to do this, however, it is necessary to understand the relationship be- 
tween the distributions of the H2 and the CO. They have rather different formation 
mechanisms: H2 forms on the surface of dust grains, while CO forms in the gas- 
phase as a product of ion-neutral chemistry. Both are readily photodissociated 
by ultraviolet (UV) radiation, but H2 can protect itself from this radiation via 
self-shielding even in relatively low column density gas (Draine & Bertoldi 1996). 
We therefore expect to find large variations in the CO/H2 ratio within any given 
GMC, with the H2 filling a significantly larger volume of the cloud than the CO. 

Despite this, there is good evidence that CO emission is a good tracer of H2 
mass within the Milky Way (see e.g. Solomon et al. 1987, Dame et al. 2001]) ■ A 
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number of independent studies have shown that GMCs in the Galactic disk show 
a good correlation between the integrated intensity of the J = 1 — > rotational 
transition line of 12 CO and the H 2 column density. This correlation is typically 
described in terms of a conversion factor Xco (the 'X- factor'), given by 

X co = ~ 2 x 10 20 cm- 2 K- 1 km" 1 s, (1.1) 

where Wco is the velocity-integrated intensity of the CO J = 1 — » emission line, 
averaged over the projected area of the GMC, and N-r 2 is the mean H2 column 
density of the GMC, averaged over the same area. 

However, the issue of the environmental dependence of Xco remains highly 
contentious. Extragalactic measurements of Xco that use a virial analysis to 
determine cloud masses find values for Xco that are similar to those obtained 
in the Milky Way, with at most a weak metallicity dependence (e.g. Rosolowsky 
et al. 2003, Bolatto et al. 2008). On the other hand, measurements that constrain 
GMC masses using techniques that do not depend on CO emission consistently 
find values for Xco that are much larger than the Galactic value and that are 
suggestive of a strong metallicity dependence (e.g. Israel H9971 Leroy et al. [2009). 

Numerical simulations provide us with one way to address this observational 
dichotomy. If we can understand the distribution of CO and H2 in realistic models 
of GMCs, then we may begin to understand why the different types of observation 
give such different results. In this contribution, we summarize the results from 
some of our recent numerical simulations that self-consistently model both the 
chemistry and the turbulent dynamics of the gas within GMCs, and discuss what 
they can tell us about the behaviour of -^co ■ 

2 Method 

We have performed a large number of simulations of the chemical and thermal evo- 
lution of the turbulent, dense interstellar medium using a modified version of the 
ZEUS-MP magnetohydrodynamical code. Our modifications include the addition 
of a simplified treatment of hydrogen, carbon and oxygen chemistry, a detailed 
atomic and molecular cooling function, and a treatment of the effects of ultravio- 
let radiation using a six-ray approximation. Full details of these modifications can 
be found in Glover et al. (|2010|) . 

Our simulations begin with initially uniform atomic gas, threaded by a uni- 
form magnetic field with strength Bq = 5.85 /iG. The initial velocity field is 
turbulent, with power concentrated on large scales, and with an initial rms veloc- 
ity of 5 kms' 1 . We drive the turbulence so as to maintain approximately the same 
rms velocity throughout the simulation. We adopt periodic boundary conditions 
for the gas and in most cases use a cubical simulation volume with a side length 
L — 20 pc. In a few simulations, we adopt a smaller-sized box, with L = 5 pc. We 
have run simulations with a variety of mean densities and metallicities, in order 
to span a range of different physical conditions. Full details of these simulations 
can be found in Glover & Mac Low (|2010|) . 
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Fig. 1. Estimate of the CO-to-H2 conversion factor Xco, est, plotted as a function of the 
mean visual extinction of the gas, (Ay). At (Ay) > 3, the values we find are consistent 
with the value of Xco — 2 x 10 20 cm~ 2 K" 1 km" 1 s determined observationally for the 
Milky Way by Dame et al. (|2001f) . indicated in the plot by the horizontal dashed line. At 
(Ay) < 3, we find evidence for a strong dependence of Xco, est on {Ay). The empirical 
fit given by Equation l3.f l is indicated as the dotted line in the Figure, and demonstrates 
that at low (Ay), the CO-to-Eb conversion factor increases roughly as Xco, eat oc Ay 3 ' 5 . 



3 Results 

Because the CO in many of our simulations is optically thick, an accurate deter- 
mination of Wco would require a full non-LTE radiative transfer calculation, a 
complex undertaking outside the scope of our present study. Instead, we make 
use of a simpler procedure to determine an estimate for the CO-to-H2 conversion 
factor, denoted as Xco, est- We first select a set of independent sightlines through 
our simulation, one per resolution element. We next compute H2 and CO col- 
umn densities along each of these sightlines. We convert each of the CO column 
densities into an estimate of the optical depth of the gas in the CO J = 1 — » 
transition, under the assumptions that (a) the CO level populations are in LTE, 
(b) the gas is isothermal, with a temperature equal to the CO-weighted mean tem- 
perature found in the actual simulation, and (c) the CO linewidth is uniform, and 
is given by Av = 3kms _1 . Given an estimate for the CO optical depth, we then 
compute an estimate for lUco using the same technique as in Pineda et al. (2008) . 
Finally, we average over all the sightlines to compute a mean intensity {W co ) for 
the simulation, and do the same for the H2 to arrive at a mean H2 column density. 
-Xco, est is then simply the ratio of these two values. 

Using this procedure, we have computed Xco. est f° r each of our simulations. 
Figure [T] shows how the values we obtain depend on the mean visual extinction of 
the gas, (Ay). We see from the figure that there is a clear change in the behaviour 
of Xco, est at (^v) ~ 3. In clouds with larger mean extinctions, Xco, est is roughly 
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constant and has a value consistent with the observationally-determined value 
for the Milky Way. On the other hand, for smaller mean extinctions, Aco.est 
increases sharply with decreasing (Ay). This behaviour is caused by a rapid 
fall-off in the CO abundance with decreasing mean extinction, which leads to 
a corresponding sharp drop in Wco ■ Because it does not self-shield efficiently, CO 
molecules are protected from photodissociation primarily by dust extinction, and 
so as this decreases, the CO abundance decreases much more rapidly than the H 2 
abundance, resulting in the rapid increase we find for Xco,cst- The dependence of 
^CO.ost on (Ay) can be described by the empirical fitting function 



illustrated in Figure Q] by the dotted line. 

The sharp fall-off in Wco with decreasing mean extinction has an important 
consequence. In order to detect CO emission from GMCs in Local Group galaxies, 
the integrated intensity of the emission must be ~ 1 K kms" 1 or higher, and we find 
in our simulations that only the clouds with (Ay) > 1 have integrated intensities 
above this value. Therefore, GMCs detectable in CO will always sit on the right- 
hand side of Figure [H in the regime where Xco.est is roughly constant, providing 
a simple explanation for why values of Xco m extragalactic systems derived using 
CO observations are always roughly the same as the Galactic value. On the other 
hand, GMCs with mass determinations that do not rely on CO are not constrained 
to fall in the regime where Wco is large, and so may occur anywhere in the plot, 
thereby explaining why they are often found to have X-factors that are much larger 
than the Galactic value. A further consequence of the behaviour of Xco and Wco 
is the prediction that CO-bright clouds in low metallicity systems must be larger 
and/or denser than their Milky Way counterparts, since at lower metallicity a 
larger surface density of gas is required to produce the necessary mean extinction. 
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